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ABSTRACT 

We present accurate models of the gravitational potential produced by a radially expo¬ 
nential disk mass distribution. The models are produced by combining three separate 
Miyamoto-Nagai disks. Such models have been used previously to model the disk of the 
Milky Way, but here we extend this framework to allow its application to disks of any 
mass, scalelength, and a wide range of thickness from infinitely thin to near spherical 
(ellipticities from 0 to 0.9). The models have the advantage of simplicity of implemen¬ 
tation, and we expect faster run speeds over a double exponential disk treatment. The 
potentials are fully analytical, and differentiable at all points. The mass distribution 
of our models deviates from the radial mass distribution of a pure exponential disk 
by <0.4 % out to 4 disk scalelengths, and <1.9 % out to 10 disk scalelengths. We 
tabulate fitting parameters which facilitate construction of exponential disks for any 
scalelength, and a wide range of disk thickness (a user-friendly, web-based interface 
is also available). Our recipe is well suited for numerical modelling of the tidal effects 
of a giant disk galaxy on star clusters or dwarf galaxies. We consider three worked 
examples; the Milky Way thin and thick disk, and a disky dwarf galaxy. 
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1 INTRODUCTION 

The mass distribution of the stellar disk of most g alaxies is 
well represented by a radially exponential profile (|Freemanl 
■ limi t. it is advantageous to be able to accurately model 
the potential field, accelerations, or tides that arise from 
a mass distribution with an exponential profile. A radially 
exponential disk profile has the following form: 

E (R) = E 0 exp(— R/Rd) (1) 

where E is the surface density, Eo is central surface density, 
R is radius within the disk, and Ra is the disk scalelength. 

Radially exponential disks may have different vertical 
density distributions. For galaxy disks, a commonly used 
form for the vertical density distribution is a sech n form: 

p(R,z) = po exp( R/Rd )sech n ( \z\j zo )• (2) 

where z o is the scaleheight, and n is typically ~1 to 3. 

* E-mail:rsmith@astro-udec.cl 


Another form of radially exponential disk, the ‘double 
exponential’, has an exponentially decaying vertical distri¬ 
bution: 

p(R, z) = po exp(-R/R d ) exp(—|«|//i z ). (3) 

where h z is the exponential disk scaleheight. In fact the 
double exponential is a special case of Eqn. [2] when 
n—>> 00 . To calculate the potential from a double exponen- 
tial disk, it is necessary t o perform the following integral 
([Binnev fe Tremainelll987h : 

$(R,z) = - —!■- / dz'exp(-z'/h z )x 

Rd J ~°° (4) 

/ oo / on/ \ V ' 

dR'srn - 1 l A+ * A ) R'K 0 (R'/R d ) 

where A + = + (R' + i?) 2 , A- = a/z 2 + (R‘ ~ 4) 2 

and Ko is the modified Bessel function of the second kind. 
This integral cannot be performed analytically, and there- 
fore is calculated nu merically (e .g iDehnen fe BinnevH l998: 
the GALPY package: iBov 3 Ena), except for in the special 
case where an infinitely thin exponential disk is assumed. 
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Due to these limitations, the potential of disk galax¬ 
ies has often b een modelled using a single Miyamoto-Nagai 
(MN) disk fe.glAllen &; Santillanlll99ll: Fellhauer et al.ll2006l : 


iFellhauer et al.ll2007l : lKupper et ahlboioUSmith et al.ll2013n . 

as this is analytical, and fully defined and provides con¬ 
tinuous derivatives at all points. The potential of a single 
Miyamoto- Nagai (MN) disk is descr ibed by the following 
expression (|Mivamoto fc Nagailll975h : 


<h(i?,z) = 


—GMmn 


U? 2 + (a + V? 2 + 6 2 ) 2 


(5) 


where Mmn is the total disk mass, a is the radial scale- 
length, and b is the vertical scaleheight. This expression can 
be trivially differentiated in the R and z direction to produce 
expressions for the acceleration at any location. If a = 0, the 
potential of a single MN disk reduces to that of a Plummer 
distribution (i.e spherical). For b = 0, the p otential reduce s 
to that of an infinitely thin Kuzmin disk (iKuzmi nl ll956l j. 
Hence by varying the parameters a and 5, mass distributions 
can be modelled with a range of thicknesses. Additionally, 
the density of a single MN disk is given by: 


p(R,z) = 


Mmn^ 2 [clR 2 + (a + 3 Vz 2 + b 2 )(a + %z 2 + b 2 ) 2 


4tt [i? 2 + (a + Vz 2 + & 2 ) 2 ] 5/2 (z 2 + 6 2 ) 3/2 


( 6 ) 

While most galaxy disks are well represented by a radi¬ 
ally exponential profile (see Eqn. [I]), a single MN disk is a 
poor match to a radially exponential disk. Its surface density 
near the centre is too low, and it attains too high surface 
density at large radius. We quantify the deviation from a 
pure exponential disk for the models considered in this pa¬ 
per in the following manner. We calculate the difference in 
mass found within a radial annuli of a single MN disk and 
a pure exponential disk, and sum up the absolute of these 
differences in all annuli out to a chosen radius. In this way, 
we calculate that the total mass deviation of a single, thin, 
MN disk from a pure exponential disk is 5.0% at 4 R d , and 
19.9% at 10 R d . 

To imp r ove th e match to a radially exponential disk, 
IFlvnn et al.l (Il996h combined three MN (3MN) disk pro¬ 
files, each with different radial scalelength a, one with neg¬ 
ative mass, and all with a single vertical scaleheight b. This 
model is better at matching a radially exponential disk 
at large radii, and we calculate the total mass deviation 
from a thin, radially exponential disk is 9.4% at 4 R d: and 
10.0% at 10 R d . This 3MN model has since been used ex¬ 
tensively for m odelling the potential from the disk of the 
Milky Way (e.g Hannin en fe Flvnnll2004l : lRodionov fe Orlovl 
12008 k iMoni Bidin et alT 2014h . However minor alterations 


to the original parameter values have been used to better 
match more recent measurements of the Milky Way ’ s disk 
scalelength and circular velocity (iGardner fe Flvrml l2010l : 
Ijflkova et alJl20l3 : iMovano Lovola et al.ll201? b 

The aim of this stu dy is to extend the utility of the 
framework introduced by I Flynn et al .1 (Il996h , a framework 
which was developed specifically for the recovery of an Milky 
Way-like radially exponential disk using 3MN potentials. 
Our extension here allows its general application to disks 
of any mass, scalelength, and thickness, rather than simply 
for the singular purpose for which it was initially designed. 
Our new models also better match the distribution of a radi¬ 


ally exponential disk. In jj2] we derive the new 3MN models, 
in jj3] we consider three worked examples, in j]4|we compare 
a 3MN model to other well known disk models, and finally 
we summarise and conclude in ^5] 


2 NEW TRIPLE MN (3MN) POTENTIALS 

From Eqn. [5] a single MN potential has 3 free parameters; 
disk mass (Ma), rad ial scaleleng th (a) , and vertical scale- 
height (5). Following IFlvnn et al (1996), we choose a single 
fixed value of b for all three MN potentials. This substan¬ 
tially reduces the total parameter space we must consider, 
and ensures our models have a uniform scaleheight. It also 
enables us to control disk thickness through a single parame¬ 
ter. Thus our 3MN models have a total of 7 free parameters; 
Mmn,i, Mmn,2, Mmn,3, ai, &2, &3, and b. Our aim is to find a 
combination of these 7 parameters that minimises the mass 
deviation from a radially exponential disk, out to 4 rR 

In order to find a good combination we used a brute 
force approach. We fix the value of 5, then numerically ran 
through a grid of values for the other 6 parameters Mmn,i, 
Mmn,2, Mmn, 3, ai, «2, a3. For each parameter set, the mass 
deviation from a pure exponential disk out to 4 and 10 R d 
was quantified numerically. Due to the significant numbers 
of combinations of parameter values possible, we initially 
ran with a coarse grid, and then later ran with a finer grid 
focussed on the best matches. 

Using this brute-force approach, we find that for an in¬ 
finitely thin disk (b/R d =0.0), a 3MN model can be found 
that radially deviates from a radially exponential disk by only 
0.4% within 4 Rd, and only 1.9% out to 10 R d . 


2.1 A recipe for different disks: varying disk 
mass, size and thickness 

Our choice of 3MN model has a highly useful property - 
once a good match to a radially exponential disk of mass 
M d and scalelength R d has been found, it can be easily 
scaled to different masses and disk sizes. For example, an 
equally good match can be found for a disk twice as massive 
by simply doubling all the MN disk masses (Mmn,i, Mmn, 2 , 
ATmn,3) • 

Therefore to calculate the best 3MN model with a differ¬ 
ent thickness, a new 6-parameter set must be found. How¬ 
ever, if the change in b/R d is very small, and smooth, we 
can expect that each parameter changes in a smooth, and 
continuous manner. We initially searched for best matches 
for different disk thicknesses using the brute force approach. 
However in practice it was difficult to find a smooth fit-line 
through these points for which every point on the fit-line pro¬ 
vided a good radial match to an exponential. This was also 
highly time consuming as the brute force approach involves 
searching through all possible combinations of a broad range 
of values for each of the 6 parameters. Thus we changed our 
approach - instead we use our best solution for an infinitely 
thin disk, found by the brute-force approach, as a prior, and 


1 The choice of 4 R d is rather arbitrary, however >90% of an 
exponential disk mass is enclosed within 4 R d 
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Figure 1. (Top panel) Evolution of the three mass parame¬ 
ters, and (middle panel) three scalelength parameters of the 3MN 
model, as a function of disk thickness b/R a - Cross symbols are the 
solution from the brute-force approach for an infinitely thin disk. 
Solid lines are the continous solutions for a range of disk thick¬ 
nesses. Dashed lines are 4th-order fits to the solid lines. (Lower 
panel) The fractional difference in mass between the model and 
a pure exponential disk measured radially outwards to 4 scale- 
lengths (black), and 10 scalelengths (red). 


Parameter k\ k 2 ks kd ks 


MyiN,l/Md 

-0.0090 

0.0640 

-0.1653 

0.1164 

1.9487 

A/MN,2/Afd 

0.0173 

-0.0903 

0.0877 

0.2029 

-1.3077 

AfMN,3/A/d 

-0.0051 

0.0287 

-0.0361 

-0.0544 

0.2242 

ai/Rd 

-0.0358 

0.2610 

-0.6987 

-0.1193 

2.0074 

a 2 /Rd 

-0.0830 

0.4992 

-0.7967 

-1.2966 

4.4441 

as/Rd 

-0.0247 

0.1718 

-0.4124 

-0.5944 

0.7333 


Table 1. Table of constants in Eqn. [TJ providing an accurate fit 
to the variation of each of the 6 parameters (see column 1) of the 
3MN model shown in Fig. \l\ 


then use an alternative approach to find the variation of each 
parameter as a function of b/Rd. 

See upper and central panel of Fig. HI The best solution 
from the brute force approach is shown by cross symbols. 
We then shift along the b/Rd axis in small steps, searching 
for new solutions that allow a continuous path across the 
figure, while simultaneously attempting to minimise the dif¬ 
ferences between the 3MN model and that of an exponential. 
In practice this approach is significantly less computation¬ 
ally challenging in comparison to the brute-force approach, 
as only a small range of possible values for each parame¬ 
ter are permitted in order to form a continuous path. As a 
result, in a fraction of the computational time required for 
the brute-force approach to find a single parameter set, hun¬ 
dreds of parameter set solutions are found for a wide range 
of b/Rd value. These solutions naturally form a continuous 
path, allowing a user to choose an arbitrary value of b/Rd 


Figure 2. Caption as in Fig [I] except now for models with pos¬ 
itive densities at all positions for the thickness range b/Rd=0 to 


1.35. 


Parameter 

ki 

k 2 

ks 

kd 

ks 

A/MN,l/Afd 

0.0036 

-0.0330 

0.1117 

-0.1335 

0.1749 

AfMN,2/Afd 

-0.0131 

0.1090 

-0.3035 

0.2921 

-5.7976 

T/mn , 3 / Ala 

-0.0048 

0.0454 

-0.1425 

0.1012 

6.7120 

ai/Rd 

-0.0158 

0.0993 

-0.2070 

-0.7089 

0.6445 

a 2 /Rd 

-0.0319 

0.1514 

-0.1279 

-0.9325 

2.6836 

as/Rd 

-0.0326 

0.1816 

-0.2943 

-0.6329 

2.3193 

Table 2. Table of constants 

in Eqn. [3 providing 

an accurate fit to 


the variation of each of the 6 parameters of the 3MN model shown 
in Fig. [2] These models have positive densities at all positions for 
the disk thickness range b/Rd from 0 to 1.35. 


solid lines in Fig. [l] and cover a range of disk thickness from 
b/Rd=0.0 to 3.0. 

The dashed lines in the upper and central panel show a 
4th-order fit to the solid lines. This fit produces 3MN models 
that match a radially exponential disk to <1.0% out to 4 Rd, 
and <3.3% out to 10 Rd (dashed black and red line in lower 
panel respectively), for the range of disk thickness considered. 
The 4th-order fit to each parameter has the form: 


Parameter — kix 4 + k 2 X 3 + ksx 2 + kdX + ks (7) 

where x = b/Rd■ The values of the constants (ki-ks) are 
given in Tab. |T] 


2.2 Issues with negative densities 

A shortcoming with the models described in Fig. [T] and 
Tab. [T] are that the MN model with negative mass, also 
has the largest scalelength. As a result it is inevitable that 
at sufficiently large radius, negative densities will be found. 
Negative de nsities at large radii also ocurred in the original 
iFlvnn et al.l (1996) model. In practice, the negative densities 


within the permitted range. The solutions are shown by the 


occur near the plane of the disk, and only in the very outer 
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Figure 3. Rotation curve of an infinitely thin exponential disk 
with Md=1.0xl0 10 M®, and Rd= 1-0 kpc (black dashed line). For 
comparison, the rotation curve of the 3MN models from Tab. 0 is 
shown (red solid line), and the rotation curve of the 3MN models 
from Tab [2] is shown (blue dash-dot line). Within 4 scalelengths, 
all models are barely distinguishable. 


disk (atR=5.6toll.2 scalelengths for b/Rd in the range 0.0 
to 3.0). As a result the negative densities are very small, and 
if the disk is placed within a dark matter halo, then the neg¬ 
ative densities will be more than offset by the halo leaving 
only positive densities everywhere. If this is of concern (e.g if 
disks without halos are studied), we also present alternative 
models. These models are advantageous as they have posi¬ 
tive densities at all positions, for disk thickness b/Rd from 0 
to 1.35. For b/Rd >1.35, negative densities begin to appear 
in the outer disk of these models as well. The lower panel of 
Fig. [2] shows that these models are still very accurate out to 
4 scalelengths (solid black line): <1% difference from a pure 
exponential). However they are often less accurate out to 
10 scalelengths (solid red line): 1.5% to 7% difference from 
a pure exponential, with the value being quite sensitive to 
disk thickness. Once again the dashed lines show a 4th or¬ 
der fit of the form shown in Eqn. 0 We provide the values 
of the constants Aq-fe for these additional models in Tab. 
[2j Although we also found other models with positive den¬ 
sities at all points for disks as thick as 5/Rd ~ 2.0, these 
models were significantly less accurate matches to radially 
exponential disks, and we will not consider them further. 


2.3 Rotation curves with the 3MN models 

In Fig.0 we compare the rotation curve for an infinitely thin 
exponential disk (black dashed line) to those of the 3MN 
models from Tab. \l\ (red solid line) and from Tab. [2] (blue 
dot-dash line). Within 4 scalelengths, the rotation curves are 
all virtually indistinguishable. At large radius the greater 
accuracy of the models in Tab. [T| is visible in comparison 
to the Tab. [2] models. We note that the negative densities 
arising in the Tab. [I] models, appear in the plane of the disk 
no closer than 8.0 scalelengths for an infinitely thin model. 
They have negligible effect on the rotation curve. 



Figure 4. Conversion from disk thickness {b/Rd) to ellipticity 
for 3MN models viewed edge-on, and with ellipticity measured 
for contours containing 25% (black), 50% (red), and 75% (green) 
of the total mass. The disk ellipticity is found to be roughly equal 
for all three contours. A quadratic fit is made to all the data points 
(blue dashed line) and the fit is given in Eqn. [8] 

2.4 Useful conversion formulae 

2.4.1 Disk thickness (b/Rd) and ellipticity (e) 

To aid comparison with observed disk ellipticities, we mea¬ 
sure how the ellipticity e of our disk models, seen edge-on, 
varies with disk thickness b/Rd. We measure the ellipticity 
of contours containing 25%, 50%, and 75% of the total mass, 
for a range of disk thickness from b/Rd = 0 to 3. The results 
are shown in Fig. 0 In practice, the ellipticity changes very 
little whether we use the contour containing 25%, 50% or 
75% of the total mass. This indicates the ellipticity of the 
3MN model is only a weak function of radius. We fit the 
data points with a quadratic formula, shown by the blue 
dashed line in Fig. 0 giving ellipticity e as a function of 
disk thickness 5/Rd: 

e = -0.099(5/R d ) 2 + 0.599(5/R d ). (8) 

Solving Eqn. 0 we obtain an expression for converting 
from disk thickness 5/Rd as a function of ellipticity e: 


5/R d = 3.025 - 3.178V0.906 - e. (9) 


2/.2 Disk thickness to exponential and sech 2 scaleheight 

Our 3MN models have a single parameter controlling their 
thickness, 5/Rd- However users are likely to be unfamiliar 
with 5/Rd as a measure of thickness. For this reason, we 
characterise this to scalelengths for more familiar vertical 
density distributions; an exponential, or a sech 2 decay with 
z. 

We first sum up the fractional difference between the 
density of the 3MN model and the exponential, measured 
vertically out of the plane up to 5 scaleheights (z =0 to 55), 
measured at R=0. We vary the disk thickness 5/Rd in order 
to minimise the sum, and best match the exponential distri¬ 
bution. We then tabulate the best matches we find between 
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Figure 5. Conversions between 3MN disk thickness (b/Rd) and 
(left panel) exponential scaleheight, or (right panel) sech 2 scale- 
height. Data points are found using the technique described in 
Sect. 12.4.21 The equations provide smooth fits through the data 
points. In each panel, the upper-left equation gives conversion to 
b/Rd, and the lower-right equation gives conversion from b/Rd- 
Conversions are approximate, but are best matches. 


b/Ra and h z /Rd over a range of b/Rd from 0 to 3. We repeat 
this procedure to find the best match between 3MN disks 
and a sech 2 vertical distribution. The best match is the same 
if we instead choose to compare the vertical density profiles 
at R = 2 or 4 Rd- 

The tabulated best matches are plotted in Fig. [5j The 
given formulae provide smooth fits through the tabulated 
data points. In each panel, the upper-left equation provides 
conversion to b/Rd, while the lower-right equation provides 
conversion from b/Rd- We emphasise that our 3MN models 
are mathematically ill-equipped to tightly match an expo¬ 
nential or sech 2 vertical density distribution at all z (an ex¬ 
ample will be shown in Sect. 0 . Thus the purpose of these 
equations is only to provide users with the approximate (but 
best available) match. This in turn provides a sense of the 
physical meaning of a particular choice of b/Rd value. 


3 WORKED EXAMPLES 

Finally we consider three worked examples; the thick disk of 
a disky dwarf galaxy, the thin disk of the Milky Way (MW), 
and the thick disk of the MW. Projected density plots for 
all three models are shown in Fig. [ 6 ] For these examples we 
use the more accurate Tab. [T] models. 


3.1 A dwarf disk galaxy with a thick disk 

We consider a dwarf disk with a total disk mass M d = 
1 x 10 8 Mo, and Rd — 0.25 k pc (iFat hi et al. l2010h . We 
assume an ellipticity of 0.6 (| Sanchez-Janssen et al.ll 2010 h . 
From Eqn. O b/Rd = 1.27. Substituting b/Rd into Eqn. 0 
with constants from Tab. [I] we get; Mmn,i = 1.94 x 10 8 M©, 
Mmn,2 — —1.05 x 10 8 Mq, Mmn,3 — 0.142 x 10 8 M©, 
a\ = 0.29 kpc, 0-2 = 0.58 kpc, CL 3 = —0.10 kpc. This model 
matches a radially exponential disk to better than 0 . 8 % at 
4 Rd, and 3.0% at 10 Rd- 


3.2 The thin disk of the MW 

We consider a radially exponential disk for the MW’s thin 
disk with a total mass of Md = 4.6xl0 10 M©, and ra¬ 
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Figure 6. Projected surface density plots for (upper) the thin 
disk model, centre the thick disk model, (lower) the dwarf galaxy 
model. The left column is face-on projection, and the right column 
is edge-on projection. The grey scale bar gives surface densities 
in units of M© pc -2 . All panels are 10 kpc on a side, to give a 
sense of relative size. 


dial exponential scalelength Rd = 2.2 kpc ([Bow fe Rixl 
§013). We choose an expone ntial scaleheight of h z = 0.2 kpc 
(lLarsen & Humphreydl2003 ). and Fig. Ogives b/Rd — 0 . 11 . 
Substituting b/Rd into Eqn. [7] with constants from Tab. [T| 
we get Mmn,i — 9.01 x 10 10 M©, A7mn,2 — —5.91 x 10 10 M©, 
Mmn,3 — 1.00 x 10 10 M©, ai = 4.27 kpc, <22 = 9.23 kpc, 
<23 = 1.43 kpc. This model matches a radially exponential 
disk to better than 0.5% at 4 Rd, and to better than 1.8% 
at 10 Rd- 


3.3 The thick disk of the MW 

For the thick disk, we consider a radially exponential disk 
with scalelength Rd — 3 .8 kpc and scaleheight h z = 0.9 kpc 
([Moni Bidin et al.ll2012h . Using Fig. Ogives b/Rd = 0.30. 
We assume that t he thick disk mass i s 8 . 6 % of the thin 
disk mass (lYoachim fe Dalcantonl 12006 ) so the total disk 
mass Md = 4.0 x 10 9 M©. Substituting b/Rd into Eqn. 0 
with constants from Tab. Ogives; M m n,i = 7.88x 10 9 M©, 
Mmn,2 — —4.97x 10 9 M©, Mmn,3 — 0.82x 10 9 M©, 
a\ — 7.30 kpc, cl 2 = 15.25 kpc, a 3 = 2.02 kpc. This model 
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Face-on surface density comparison 



the plane (right). The scaleheight of the exponential and 
sech 2 distribution is the best match to the 3MN model, cal¬ 
culated using the equations in Fig. [5] Clearly it is impossible 
for the 3MN model to exactly match the other profiles, as 
they are mathematically distinct. However, by comparing 
them we can judge how they differ, and the quality of the 
best matches provided in Fig. [5] Comparing along each row, 
the exponential disk is most ‘cuspy’ in the vertical direction, 
where as the sech 2 is the most ‘cored’, and the 3MN model 
is found somewhere in between the other two. This is true 
whether the cross-section is made at (a) y=0Ra , (b) y=2Rd , 
or (c) y=AR d . 

Row (d) provides a more quantitative analysis of the 
vertical density distribution up to ~56, measured at (left) 
R=0Rd, (centre) R=2Rd, and (right) R=4Rd- At small z 
(less than roughly 2 or 36) the 3MN density distribution 
roughly matches the other two profiles. However at greater 
distances from the plane, the 3MN model returns higher 
densities, and this becomes stronger if the vertical density 
profile is measured at larger R. 

The upper panel of row (e) compares the surface density 
of the 3MN model (solid curve) to a radially exponential disk 
(dashed curve). Curves are plotted out to R=4Rd- There is 
clearly an excellent agreement over this radius range, al¬ 
though it can be seen that the 3MN model returns slightly 
lower densities at very small radii (R<0.2Rd). The lower 
panel indicates the absolute fractional difference in surface 
density between the two profiles. The maximum density dif¬ 
ference is ~15% at R— 0, but for R>0.2Rd the density dif¬ 
ference is very low (typically ~0.5%). 


5 SUMMARY AND CONCLUSIONS 


Figure 7. Analysis of the MW thick disk model of Sect. 13.31 
Rows (a), (b), and (c) are cross-sections through the volume den¬ 
sity distribution, and each row allows comparison between the 
3MN model (left), a double exponential (middle), and a radially 
exponential model with a sech 2 density drop off out of the plane 
(right). Grey scale color bar units are Mq pc~ 3 . The plane of 
the disk lies on the x-y plane, and slices are perpendicular to the 
plane at (a) y=0Rd, (b) y=2R d, and (c) y=AR^- Row (d) is vol¬ 
ume density with distance from the disk plane, measured at (left) 
ORd, (centre) 2Rd, an d (right) AR&- Different curves (indicated 
in key) compare the 3MN model, double exponential, and sech 2 
model. The upper panel of row (e) compares the surface density 
profile of the 3MN model (solid curve) with a radially exponential 
disk (dashed curve) out to AR&. The lower panel of row (e) shows 
the absolute fractional difference in the surface density of the two 
profiles at each radius. 

matches a radially exponential disk to better than 0.6% at 
4 Rd, and to better than 1.4% at 10 Rd- 


4 COMPARISON OF 3MN MODEL WITH 
OTHER DENSITY PROFILES 

In Fig. [71 we analyse the 3MN thick disk model from Sect. 
f373l in more detail. In row (a)-(c) we compare cross-sections 
through the volume density distribution of the 3MN model 
(left), a double exponential (middle), and a radially expo¬ 
nential model whose density profile decays as sech 2 out of 


We present a recipe for using a triple Miyamoto-Nagai 
(3MN) disk distribution to model the potential of a radially 
exponential disk. 3MN models have previously been used to 
model the disk component of the Milky Way. Here we ex¬ 
tend on this framework to allow its general application to 
disks of any mass, scalelength, and a wide range of thick¬ 
ness. We find parameters of the 3MN model that best match 
the mass distribution of a radially exponential disk. We con¬ 
sider a broad range of disk thicknesses from infinitely thin 
(ellipticity = 0.0), to near spherical (ellipticity = 0.9). The 
3MN models have many benefits as they are entirely analyt¬ 
ical, easy to implement, and provide continuous derivatives 
(enabling a calculation of accelerations) at all points. 

(i) We provide accurate fitting formulae to our new 3MN 
models, that reproduce the mass distribution of a radially 
exponential disk to <1.0% out to 4Rd, and <3.3% out to 
10Rd for disks with a range of ellipticities from flat to near 
spherical (see Eqn.[7]and Tab. [l]). 

(ii) We provide a second set of models in Tab. [2] that en¬ 
sures positive densities at all positions for the disk thickness 
range b/Rd = 0.00 to 1.35. This is equivalent to an ellipticity 
range from 0.0 to 0.6. 

(iii) We provide a fitting formula to allow for easy con¬ 
version between the disk thickness 6/Rd and disk ellipticity 
e (see Eqn.[8]), and in reverse (see Eqn.(9]). 

(iv) The vertical distribution of our 3MN models is sim¬ 
ilar to an exponential or sech 2 distribution at small z. We 
provide a rough approximation for converting between disk 


















thickness, exponential scaleheight, and sech 2 scaleheight for 
z up to 5 scaleheights (see equations given in Fig. [5]). 

A user-friendly, online web-form is available at 
http://astronomy.swin.edu.au/~ cflynn/expmaker.php. 
Users can request the disk thickness they require. The 
page will automatically provide the best matching 3MN 
parameters, calculated using our scheme. We acknowledge 
that better 3MN solutions may exist, hidden in the param¬ 
eter space. In the future we will extend on the techniques 
developed in this study in order to uncover alternative and 
potentially more accurate 3MN parameter sets, but also to 
allow for flared disks, and for alternative vertical density 
distributions. 
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